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INTRODUCTION 


The  surface  morphology  of  a  polyurethane  coating  is  an  important  factor  in  determining  a 
coating's  appearance  and  properties.  Extended  UV  exposure  can  alter  this  morphology  by 
breaking  down  the  resin  that  binds  the  fine  paint  particulates.  Figure  1  shows  scanning  electron 
microscope  (SEM)  images  illustrating  the  effect  of  UV  exposure  on  a  polyurethane  coating.  The 
figure  shows  the  resin  degradation  that  changes  the  appearance  of  the  polyurethane  coating  under 
ambient  light.  There  are  two  major  effects  of  resin  on  appearance.  First,  the  relative  amount  of 
resin  in  the  paint  controls  the  degree  of  gloss.  Second,  the  resin  serves  to  promote  a  uniform 
index  of  refraction  throughout  the  filler  particle/resin  system  in  order  to  provide  sufficient 
transparency  so  that  the  dye  particles  are  observable.  Resin  degradation  can  thus  affect 
appearance  by  altering  the  surface  gloss  through  surface  roughening  and  by  promoting  scatter 
from  the  loosened  particles.  Both  factors  can  contribute  to  the  UV  fading  (lightening)  effect 
observed  under  ambient  light  for  the  solvent-based  polyurethane  employed  in  the  present  study. 
However,  a  UV-exposed  pol5airethane  coating  appears  darker  than  an  unexposed  coating  when 
both  surfaces  are  coated  with  a  conductive  metal. 


Unexposed  Exposed 

Figure  1.  Scanning  electron  microscope  images  of  a  polyurethane  coating 
before  and  after  18  weeks  of  accelerated  UV  exposure. 

Laser  scanning  confocal  microscopy  (LSCM)  and  atomic  force  microscopy  (AFM)  have 
been  used  for  measures  of  specific  paint  surface  features  (refs  1,2).  In  the  present  study,  LSCM 
and  AFM  were  applied  to  quantify  surface  structures  in  terms  of  scaling  properties.  More 
specifically,  the  coating  structure  of  the  polyurethane  paint  surface  was  analyzed  to  establish 
correlations  between  the  intrinsic  spatial  scaling  properties  of  the  surface  and  UV  exposure  time. 
The  surface  was  analyzed  as  self-affine  structure  with  statistical  invariance  under  anisotropic 
dilations  from  the  10-|xm  AFM  to  the  100-p.m  LSCM  scan  sizes  employed  in  the  measurements. 
The  scaling  exponents  that  characterize  a  self-affine  surface  (ref  3)  were  initially  determined 
using  the  height-height  correlation  function.  However,  inhomogeneities  in  the  surface  structure, 
variations  in  the  curvature  of  the  sample,  and  contributions  from  dominant  features  of  the 
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substrate  have  a  significant  effect  on  this  type  of  analysis.  The  effect  is  most  noticeable  in  the 
analysis  of  the  AFM  images,  where  the  curvature  of  the  particles  introduces  inherent  bias  in  the 
data.  The  influence  of  these  effects  can  be  minimized  by  analyzing  the  surface  relative  to  a 
curved-reference  surface  (refs  4,5)  or  dividing  the  surface  into  patches  that  are  analyzed 
separately  using  the  "min-max"  method  (ref  6).  Meisel  et  al.  (ref  4)  suggest  a  technique  for 
measuring  scaling  properties  relative  to  a  curved-reference  surface  by  limiting  analysis  to  local, 
statistically  homogenous  regions  that  are  large  enough  to  define  a  statistically  reliable  height- 
height  correlation  function.  However,  determining  the  size  of  the  regions  that  properly  adapt  to 
the  heterogeneity  of  the  surface  requires  subjective  investigator  input.  The  "min-max"  method 
requires  the  investigator  to  judge  the  number  of  local  maxima  and  minima  of  the  surface  to  use 
in  the  analysis.  A  DFA  does  not  require  subjective  input  and  is  a  systematic  procedure  that  is 
well  suited  for  measuring  the  subtle  changes  in  scaling  exponents  of  paint  surfaces  that  result 
from  UV  exposure. 

Peng  et  al.  (ref  7)  introduced  DFA  to  determine  long-term  correlations  embedded  in  time- 
series  data.  The  sequence  of  integrated  data  of  length  N  is  divided  into  NA  nonoverlapping 
segments  of  length  /.  The  average  local  variance  of  the  segments  is  measured  relative  to  a  linear 
least-squares  fit  to  the  data  within  each  segment  /.  The  dependence  of  the  variance  on  I  defines 
the  scaling  properties  of  the  time-series. 

APPROACH 

In  the  present  study,  DFA  was  used  to  determine  the  scaling  properties  of  surfaces  in 
three-space  using  spatial  averages  of  roughness  for  surface  patches  of  width  I  measured  relative 
to  a  plane-fit  (Appendix  A)  to  the  local  data.  The  average  local  roughness,  is  given  by 

M'!(/,l)  =  ([A(JC.O-A/(x,r)P)^  (1) 

where  h{x,t)  is  the  height  of  a  single-valued  surface  at  location  x  relative  to  a  plane-fit  to  the 
local  data  at  time  t,  h^x,t)  is  the  average  height  of  the  local  data  at  time  t,  t  is  the  exposure  time, 
and  I  is  the  size  of  local  surface  patches.  WL(l,t)  is  measured  for  geometrically  increasing  values 
of  /(/  <  system  size  L)  for  each  exposure  time  t.  The  dependence  of  the  local  roughness  on  / 
gives  rise  to  the  scaling  properties  of  the  fluctuations  characterized  by  the  roughness  exponent  cx 
at  time  t  as  (ref  8) 


(2) 


Figure  2  shows  filtered  LSCM  data  and  the  patchwork  of  local  reference  surfaces 
employed  by  DFA  for  /  =  1.25, 2.5,  and  5.0-|j,m.  The  system  size  (scan  size)  is  L  =  10-|i.m  for 
the  AFM  data  and  L  =  100-|i.m  for  the  LSCM  data. 
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/  =  2.5^m  l=5.0]im 


Figure  2.  DFA  reference  patches  for  a  lO-iim  filtered  LSCM  image 
of  a  solvent-based  polyurethane.  (All  dimensions  in  jim.) 

The  paint  in  the  present  study  is  a  solvent-based  polyurethane  system  (MIL-PRF-85285) 
using  nonpolymeric  flattening  agents  (silica),  and  is  pigmented  to  Air  Force  gray  color  number 
37375.  The  accelerated  coating  degradation  was  conducted  using  a  QUV  weathering  chamber 
according  to  test  method  ASTM  G  53-96.  A  UVA  340  light  source  was  used  that  emitted  a 
spectral  irradiance  of  0.77  W/m^  measured  at  340-nm. 

The  time  evolution  for  each  of  the  exposed  paint  surfaces  is  not  available,  however, 
samples  are  available  for  three  different  UV  exposure  times.  The  DFA  for  these  samples 
indicates  that  wi(l,t)  does  not  saturate  at  the  maximum  exposure  time.  Therefore,  the  time 
evolution  of  the  intrinsic  surface  structure  as  a  function  of  UV  exposure  time  was  measured  in 
terms  of  the  growth  exponent,  (ref  6).  fi  characterizes  the  time-dependent  dynamics  of  the 
roughening  process  for  times  much  less  than  a  saturation  time,  as 

v’dlt)  ~  1^  (t«td  (3) 

Both  LSCM  and  AFM  were  employed  to  map  surface  structure  on  a  scale  from  100-nm 
to  100-|i.m.  These  techniques  generate  surfaces  in  three-space  that  quantify  surface  morphology 
with  little  or  no  sample  preparation.  The  AFM  is  a  type  of  scanning  probe  microscopy  that 
provides  high-resolution  images  by  scanning  a  sharp  probe  (2-nm  to  20-nm  tip  radius)  mounted 
on  a  flexible  cantilever  across  a  surface.  A  laser  is  focused  on  the  tip  and  the  beam  is  reflected 
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back  to  a  photodiode.  Deflections  in  the  cantilever  are  measured  at  the  photodiode  and  recorded 
as  the  tip  scans  across  the  surface.  The  AFM  images  in  this  investigation  have  a  20-nm 
horizontal  resolution  and  a  0.1 -nm  vertical  resolution. 

The  LSCM  is  an  optical  technique  that  generates  a  surface  three-space  using  a  series  of 
two-dimensional  images  collected  sequentially.  Each  two-dimensional  image  represents  the 
measured  reflected  intensity  of  a  laser-scanning  beam  focused  on  a  matrix  of  points  in  a  single 
focal  plane.  The  focal  plane  is  changed  incrementally  and  the  height  of  the  surface  at  each  point 
in  the  matrix  is  given  by  the  maximum  value  of  each  two-dimensional  image  at  that  point.  The 
wavelength  of  light  and  the  numerical  aperture  of  the  lens  limit  the  lateral  resolution.  The 
LSCM  images  in  this  investigation  have  a  0.25-|J,m  lateral  resolution  and  a  0.1-|i.m  vertical 
resolution  (ref  9). 

An  inherent  limitation  of  practical  LSCM  is  the  contamination  of  images  with  noise 
spikes.  These  spikes  are  a  result  of  low  surface  reflectivity,  the  objective  lens,  and  edge  effects 
of  the  surface.  They  are  a  consequence  of  associating  noise  on  a  low-level  signal  with  the 
maximum  reflected  intensity  at  a  given  focal  plane.  The  height  at  a  given  location  directly 
correlates  with  the  maximum  intensity  of  the  reflected  laser  detected  by  the  photodiode,  which  is 
assumed  to  be  at  the  appropriate  focal  plane.  However,  the  intensity  of  the  reflected  beam  may 
remain  low  at  a  given  point  in  the  two-dimensional  matrix  for  the  entire  series  of  focal  planes  in 
the  scan  (usually  256).  In  this  case,  the  height  recorded  at  that  point  may  correspond  to  a  noise 
signal  that  is  a  maximum  at  any  of  the  focal  planes  during  the  scan.  The  low  intensity  of  the 
reflected  beam  can  result  from  the  inability  of  the  scanning  beam  to  properly  focus  on  edges  of 
surface  inhomogeneities.  This  effect  introduces  spurious  spikes  around  such  edges.  Our  results 
suggest  that  as  the  complexity  of  a  surface  increases,  there  is  a  proportional  increase  in  the 
density  of  this  spike  noise. 

RESULTS 

Inspection  of  these  specimens  under  ambient  light  shows  the  familiar  feature  that 
alterations  in  surface,  appearance  vary  directly  with  UV  exposure  time.  Quantitative 
measurements  of  reflectivity  properties  are  in  progress.  The  as-deposited  paint  color  is  medium 
gloss,  light  gray,  which  becomes  significantly  whiter  in  color  under  ambient  light  with  UV 
exposure  time.  By  contrast,  the  specimens  sputter  coated  with  palladium  show  the  reverse  effect 
in  that  the  UV  exposed  specimens  darken  under  ambient  light  relative  to  an  unexposed 
palladium-coated  paint  specimen.  The  following  results  relate  directly  to  the  origin  of  these 
appearance  changes.  All  quantitative  measures  represent  averages  of  five  different  scans  for 
exposed  surfaces  and  fifteen  different  scans  for  unexposed  surfaces. 

Figure  3  is  a  comparison  of  10-p.m  AFM  and  LSCM  images  (filtered  and  unfiltered)  of  a 
typical  paint  surface.  The  figure  also  shows  the  results  of  filtering  the  LSCM  data  with  a 
nonlinear  median  filter  using  different  size  convolution  kernels.  The  median  filter  is  well  suited 
for  removing  the  spike  noise  associated  with  LSCM  images,  since  edge  information  is  preserved. 
The  paint  surfaces  shown  in  Figure  3  were  coated  with  60-nm  of  sputtered  palladium  to  prevent 
the  I^CM  scanning  laser  from  penetrating  the  paint  surface  and  ensure  the  images  represent  the 
same  surface  topology  as  measured  by  the  AFM.  The  sputtered-palladium  process  is  expected  to 
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accurately  replicate  the  substrate  for  these  substrates.  A  DFA  study  of  AFM  images  has  shown 
that  60-nm  of  sputtered  palladium  has  no  detectable  effect  on  the  surface  morphology.  Figure  4 
shows  a  comparison  of  the  coated  and  uncoated  LSCM  DFA  results  for  exposed  and  unexposed 
surfaces  using  5x5  median  filtering.  The  results  show  the  optical  surface  (uncoated)  generated 
by  the  LSCM  relative  to  the  actual  physical  surface  (coated).  The  maximum  difference  between 
the  physical  surface  and  the  optical  surface  is  near  the  O.l-pm  LSCM  resolution. 


0  0 


3x3  median  filtered  LSCM  image  5x5  median  filtered  LSCM  image 

Figure  3.  Baseline  10-|xm  AFM  and  LSCM  raw  and  filtered  images  of  a  coated, 
unexposed  solvent-based  polyurethane.  (All  dimensions  are  in  pm.) 
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Figure  4.  DFA  results  for  LSCM  images  of  uncoated  and  coated  paint 
samples  using  5x5  median  filtering.  (All  dimensions  are  in  pn.) 

In  addition  to  providing  a  means  for  measuring  the  scaling  properties  of  an 
inhomogeneous  surface  structure,  DFA  was  also  used  to  determine  the  appropriate  technique  for 
minimizing  the  considerable  effects  of  the  noise  associated  with  LSCM  images.  Figure  5  shows 
the  results  of  DFA  applied  to  AFM  and  LSCM  data  sets  collected  from  a  palladium-coated 
sample  with  18  weefe  of  UV  exposure.  A  DFA  was  employed  to  determine  the  optimal  filter  for 
removing  surface  artifacts  from  LSCM  images.  The  DFA  results  are  shown  for  raw  AFM  and 
LSCM  data  and  for  the  LSCM  data  after  median  filtering  with  different  size  convolution  kernels. 
The  lO-pm  images  shown  in  Figure  3  suggest  that  LSCM  data  that  have  been  filtered  with  a 
nonlinear  median  filter  are  consistent  witfi  the  baseline  AFM  image  considering  the  lower 
resolution  of  the  LSCM  data.  The  DFA  results  given  in  Figure  5  show  that  median  filtering  of 
LSCM  data  with  a  5x5  or  7x7  convolution  kernel  is  the  best  compromise  in  ensuring  that  the 
LSCM  data  are  consistent  with  AFM  results.  To  minimize  the  effect  of  offset  in  wl  that  can 
occur  while  acquiring  AFM  and  LSCM  images,  the  difference  in  DFA  results  {{wdiltf  - 
)  between  18- week  and  unexposed  paint  surfaces  was  also  computed  to  aid  in 
selecting  the  appropriate  filter.  Figure  6  shows  these  results  for  different  size  convolution 
kernels.  The  results  suggest  that  median  filtering  of  LSCM  images  with  a  5x5  convolution 
kernel  produces  results  consistent  with  the  baseline  AFM  data.  Therefore,  a  5x5  median  filter 
was  applied  to  all  of  the  lO-pm  LSCM  images  in  this  study. 
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Figure  5.  Local  scaling  for  coated  samples  after  18  weeks  of  exposure.  Results 
are  shown  for  unfiltered  AFM  images  and  for  filtered  LSCM  images  with 
different  size  convolution  kernels.  (All  dimensions  are  in  |im.) 
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Figure  6.  Difference  in  local  scaling  between  coated  samples  after  18  weeks  of  exposure  and 
coated  unexposed  samples.  Results  are  shown  for  the  unfiltered  AFM  images  and  for  filtered 
LSCM  images  with  different  size  convolution  kernels.  (All  dimensions  are  in  |xm.) 
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Figure  7  shows  the  results  of  DFA  applied  to  filtered  LSCM  images  of  unexposed 
surfaces  and  for  surfaces  after  3, 6,  and  18  weeks  of  UV  exposure.  The  effects  of  UV  exposure 
are  greatest  at  small  /(<  5-jim)  and  have  no  discemable  effect  on  the  surface  morphology  for  I  > 
10-|a,m.  The  results  also  show  that  a  single  exponent  cannot  adequately  describe  the  scaling 
behavior  of  the  surfaces.  Therefore,  we  employed  a  new  procedure  to  determine  if  this  type  of 
coating  structure  is  consistent  with  bimodal  behavior  having  two  unique  scaling  regions,  and  to 
obtain  accurate  measures  of  the  two  regions.  We  developed  this  technique  to  introduce  a 
systematic  method  for  extracting  the  scaling  parameters,  since  there  is  no  universally  agreed 
upon  procedure  for  estimating  the  exponent  values,  and  blind  regression  fits  to  the  height 
correlation  data  result  in  misleading  values  in  the  exponent  (ref  6).  This  technique  determines 
the  scaling  parameters  by  fitting  two  linear  and  one  polynomial  spline  to  the  DFA  results  and 
estimating  the  scaling  measures  using  the  data  that  minimize  the  residuals  in  the  fit  (ref  10).  The 
procedure  is  summarized  in  Appendix  B.  The  bimodal  fit  to  the  data  resulted  in  scaling  regions 
between  approximately  1.5  to  6-|im  and  between  5  and  15-|xm.  Analysis  was  limited  to  the  1.5 
to  5-|j,m,  where  WL(l,t)  increases  with  UV  exposure  time.  Beyond  6-|j.m,  DFA  results  overlap 
and  no  significant  effects  of  UV  exposure  are  evident. 


Figure  7.  Local  scaling  of  uncoated,  filtered,  10-pm  LSCM  images  of  a  solvent-based 
polyurethane  with  increasing  UV  exposure.  (All  dimensions  are  in  |Lim.) 
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Table  1  summarizes  the  results  of  the  bimodal  analysis.  The  roughness  of  the  surface  on 
the  1.5  to  6-|im  scale  increases  (or decreases)  monotonically  with  increased  exposure  time. 
Figure  7  shows  that  w\Xl,t)  does  not  saturate  with  time,  therefore,  the  growth  exponent,  fi,  was 
measured  by  analyzing  the  time  evolution  of  the  surface  assuming  times  shorter  than  the 
saturation  crossover  time,  tx..  P  was  determined  using  unfiltered  data  by  assuming  that  spike 
noise  is  proportional  to  surface  roughness.  The  time  evolution  wi{l,{)  was  measured  relative  to 
(variances  of  uncorrelated  data  are  additive),  which  is  nonzero  for  the  unexposed  surface. 
P  was  determined  for  values  of  I  up  to  the  limit  of  the  first  scaling  region  (6-|iim)  and  the  results 
averaged,  assuming  a  fixed  a.  This  ensured  the  evolution  of  roughness  was  a  function  of  UV 
exposure  time  only.  The  time  evolution  of  the  surface  described  by  P  was  determined  to  be  0.25 
±0.02. 


Table  1.  Bimodal  Analysis 


Exposure  Time 
(Weeks) 

First  Scaling  Region 
(|im) 

Second  Scaling  Region 
(Itm) 

a 

0 

1.5  to  5 

5  to  15 

0.60 

3 

1.5  to  6 

6  to  9 

0.58 

6 

1.5  to  5 

5  to  13 

0.56 

18 

1.5  to  6 

5  to  15 

0.53 

DISCUSSION  and  CONCLUSIONS 

The  results  in  Figures  4  through  7  show  that  there  is  generally  a  spectmm  of  local 
roughness  values  to  consider  when  quantifying  the  effect  of  UV  exposure  on  the  structure  of  the 
paint.  The  surface  structure  of  the  solvent-based  polyurethane  paint  in  this  study  changes  with 
UV  exposure,  but  only  by  small  amounts  (<  0.15-|J,m)  and  at  very  small  scales  (<  6-p.m).  The 
appearance  of  the  polyurethane  coating  under  ambient  light  changes  (i.e.,  becomes  lighter)  due  to 
a  breakdown  in  the  resin  that  binds  the  fine  particulates.  The  absence  of  flattening  additives  in 
the  paint  suggests  that  the  appearance  changes  originate  primarily  from  surface  roughness 
changes.  The  corresponding  results  for  palladium-coated  specimens,  where  light  penetration  is 
precluded,  show  that  light  scattering  associated  with  UV-induced  roughness  changes  produces 
significant  darkening.  This  illustrates  that  both  subsurface  effects  and  surface  roughness 
contribute  to  appearance  changes  associated  with  UV  exposure. 

The  resolution  of  LSCM  images  is  adequate  for  measuring  general  changes  in  surface 
morphology  of  the  paints;  however,  inherent  noise  at  scales  in  the  range  of  interest  must  be 
removed.  Using  AFM  provides  essential  reference  data  for  determining  the  best  technique  for 
removing  this  noise.  A  single  growth  exponent,  p,  can  describe  the  time-dependent  dynamics  of 
the  UV-induced  roughening  process  assuming  LSCM  spike  noise  varies  directly  with  surface 
roughness. 
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Detrended  fluctuation  analysis  is  a  valuable,  systematic  technique  for  obtaining  the 
scaling  exponents  of  an  inhomogeneous  structure.  It  is  particularly  convenient  when  analyzing 
real  data  for  scaling  properties  where  the  underlying  background  structure  is  a  dominant  factor  in 
the  measurements,  hi  the  present  study,  DFA  provided  a  means  for  quantifying  the  subtle 
changes  that  occur  in  surface  morphology  as  the  resin  deteriorates  under  UV  exposure.  The 
intrinsic  structure  of  the  paint  in  this  study  is  consistent  with  a  self-affine  surface  fractal  having 
two  distinct  scaling  regions.  Only  the  small  scaling  region  is  affected  by  UV  exposure,  and 
increases  in  complexity  with  time. 

SUMMARY 

Changes  in  the  intrinsic  structure  of  paint  surfaces  resulting  from  extended  UV  exposure 
can  significantly  alter  the  appearance  of  the  paint  due  to  a  breakdown  in  the  resin  that  binds  the 
fine  paint  particulates.  These  changes  were  characterized  in  terms  of  scaling  exponents  by 
quantifying  the  local  roughness  using  DFA  to  identify  long-range  power-law  correlations  and 
correct  for  inhomogeneities  in  the  surface  structure.  The  DFA  was  shown  to  provide  a  valuable, 
systematic  technique  for  obtaining  the  scaling  exponents  of  an  inhomogeneous  structure.  It  was 
particularly  convenient  for  the  paints  in  this  study  where  the  underlying  background  structure 
was  a  dominant  factor  in  the  measurements.  The  complexity  of  the  surface  was  shown  to 
increase  monotonically  with  UV  exposure  time,  and  the  dynamics  of  this  roughening  process 
was  described  by  a  single  scaling  exponent. 
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APPENDIX  A 


The  following  series  of  equations  provides  the  coefficients  of  a  plane-fit  (f(x,y)  =  a  +  bx 
+  cy)  to  the  data  z,  that  minimize  the  square  of  the  residuals  (^  =  a  +  bxi  +  cy,  -  Zif  for  all 
points,  m,  containing  a  local  surface  patch  of  length  /: 


m  m  m 

—  =  0=>a*m  +  b'^x.  +cjy,. 

111 

(Al) 

m  m  m  m 

^  0  =>  xf  +  c J  X,.  y,.  =  J  X-Z,. 

1111 

(A2) 

m  m  m  m 

— =0=>  a  J  y,.  +*Jx,.  y.-  +  cj  yf  =  J  y^z,- 
11  11 

(A3) 
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APPENDIX  B 


A  new  systematic  procedure  was  developed  to  determine  if  a  distribution  exhibits 
bimodal  behavior  and  to  obtain  accurate  measures  of  the  two  scaling  regions.  The  technique 
determines  the  scaling  parameters  by  fitting  a  two-linear  and  one-polynomial  spline  to  the  data 
and  estimating  the  scaling  measures  using  the  data  that  minimize  the  residuals  in  the  fit.  The 
first  linear  range  is  determined  by  fitting  a  straight  line./y,  to  the  first  ki  points,  and  the  second 
linear  region  is  determined  by  fitting  a  straight  line,/2,  to  the  second  k2  points.  A  polynomial, /i, 
is  then  fit  to  points  k2  through  m,  where  m  is  the  last  point  in  the  fit  offs,  kj  is  the  index  of  the 
data  point  at  which  the  two  linear  regions  meet,  and  k2  is  the  index  of  the  data  point  at  which /2 
and join.  Both  ki  and  ^2  are  variable  knots  (indices)  of  a  linear/linear  and  a  linear/curved 
spline  and  are  nonlinear  parameters  determined  by  exhaustion.  The  coefficients  of  the 
polynomials  are  evaluated  with  continuity  enforced  at  kj  and  kz,  and  the  residuals  (^)  are 
computed  over  the  entire  range.  The  values  of  ki  and  k2  for  which  the  residuals  are  a  minimum 
are  selected  as  the  indices  for  the  knots.  The  shape  of  the  height  correlation  data  suggests  that  a 
quadratic  fit  is  sufficient  to  determine  k2,  although  the  procedure  is  similar  if  higher  order 
polynomials  for  fs  are  used.  Continuity  is  enforced  for  both  knots  with  additional  d  continuity 
enforced  for  the  quadratic  fit  at  ^2.  If  either  f  or  fz  are  determined  to  have  fewer  than  three 
points,  then  it  is  assumed  that  only  one  linear  region  exists  and  the  structure  is  not  bimodal.  If 
both/;  and/2  have  fewer  than  three  points,  then  no  linear  region  exists.  If  the  structure  is 
determined  to  be  bimodal,  k;  is  assumed  to  be  the  limit  of  the  first  scaling  range,  kz  is  determined 
from  the  intersection  of  fz  and  the  saturation  limit  where  no  correlation  is  observed.  The 
procedure  is  outlined  below  as 


fj(x)  =  a+b(x-k,)  (Bl) 

f2ix)  =  c+d{x-ki)  (B2) 

f2ix)  =  e  +  f(x-k2)+gix-k2)^  (B3) 

Assuming  continuity  at  knots  kj  and  kz ,  with  additional  d  continuity  at  kz, 

f,(k,)  =  f2(kj),  fz(k,)  =  f,(k,),f,(k,)  =  fl(k,)  (B4) 


and  the  residuals  (£;)  become 


£i  =  a  +  bpii  -  yi ,  for  i  <  ki 
£i  =  a-¥  dpu  -  yi,  for  kj  <  i  <  kz 


Si  =  a  +  d(K+piz)  +  g(qiz)  -  yi ,  for  kzK  i  <  m 
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where 


Pii  =  Xi-ki,  qu  =  (xi-kif,  pi2  =  Xi-kz,,  qa  =  (xi-kzf  (B5) 

and  yi  is  the  measured  value  at  i. 

The  coefficients  of  the  polynomials  are  determined  by  minimizing  the  sum  of  the  squares 
of  the  residuals  S 


m  kj 


I  Jt,+; 


95  A  ”1  m  m 

■—  =  0^a*m-^bY,Pii+dY,Plj-\-d  ^gY^q.^=Yy. 

da  I  *,+7  [i^^i  J  I 


JiC  k,! 

—  =  0=>a'^  A;  +^E  =  E  y^Pn 

da  111 


^  =  0=>a  i;p,-7+Ep.-2+^^+  I^q.„  +  2KY{p,,  +  q,,)+K^^  + 


^  Z?i2(^+Pf2)l  =2]>'iA7  +  S3'/(^ 
^*2+/  J  A:;+y  *2+/ 


^  o  Tti  tn  \  tn  ni 

—  =  0=>a2^i2+^  2]^.7(^+P72)|  ,+  g'^qfz  ='Zyi^i2 

Og  k2+l  \_k2+l  J  k2+l  k2+l 


(BIO) 


where  K  =  x(k2)-x(ki),  and  A  =  m-kz. 
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